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The optical system used by Kandpal, 
Vaishya, and Joshi in their experiments 
on coherence-induced spectral shifts is 
analyzed theoretically. An approximate 
form for ihe cross-spectral density in 
the secondary source plane is obtained, 
and it is shown that, contrary to the as- 
sertions of Kandpal, Vaishya, and Joshi, 
Ihe corresponding complex degree of 
spectral coherence in this plane is 
wavelength dependent. After making 
some assumptions about the behavior 
of the interference filter used in the 
system, an approximate form for the 



spectrum of the light on-axis in the ob- 
servation plane is obtained. It is shown 
that the peak wavelengths of (his spec- 
trum do not agree with those reported 
by Kandpal, Vaishya, and Joshi. Possi- 
ble reasons for this disagreement arc 
discussed. 
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1. Introduction 

In 1989 Kandpal, Vaishya t and Joshi of the 
National Physical Laboratory (India) published re- 
sults of experiments in which they observed spec- 
tral shifts caused by a simple optical system [1], In 
their experiments light from the exit aperture of an 
integrating sphere was imaged by a two lens system 
which had an interference filter between the 
lenses. They found that when a small aperture was 
placed in the image plane, the peak wavelength of 
the spectrum of the light measured on-axis in the 
far zone of the aperture was shifted away from the 
peak wavelength which occurred when no aperture 
was used. 

They interpreted these shifts as being "Wolf 
shifts," i.e., coherence-induced shifts of the type 
predicted by Wolf [2-3]. Their explanation [4] for 
the occurrence of the shifts was based on their as- 
sertion (which was neither proven directly experi- 
mentally nor justified theoretically) that the 
filter-lens combination eliminated the wavelength 



dependence of the complex degree of spectral co- 
herence of the light in the image plane. If this as- 
sertion is true, then a coherence-induced change of 
the spectrum is to be expected when the light prop- 
agates from the image plane to the far zone, as was 
shown experimentally by Morris and Faklis [5] and 
Faklis and Morris [6]. Kandpal et al. [4] then 
argued that the introduction of a circular aperture 
in the image plane helped in modifying the spec- 
trum in the far zone. 

In this paper, we will analyze theoretically the 
optical system used in Ref. [1]. Many details of the 
system which did not appear in Ref. [1] were pro- 
vided to us by the group at NPL [7]. The basic 
outline of our paper is as follows. In Sec. 2 the 
optical system is described, and the basic assump- 
tions to be made throughout the paper are stated. 
In Sec. 3 an approximate form for the cross-spec- 
tral density of the light in the image plane is ob- 
tained, and it is shown that the corresponding 
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complex degree of spectral coherence is not wave- 
length independent. In Sec* 4 the spectrum of the 
light on-axis in the observation plane is investi- 
gated, and an approximate form for it is obtained, 
both for the case in which no aperture is used, and 
for the case in which the small aperture is used. It 
is shown that, for each interference filter used, the 
peak wavelength in the latter case is shifted with re- 
spect to the peak wavelength in the former case. 
However, these shifts do not agree with those ob- 
served by the group at NPL. The shifts predicted by 
our analysis are much smaller, so small as to be un- 
observable to within the accuracy of their experi- 
ments. In Sec- 5 our conclusions are presented. 

2. Optical System, Notation and Assump- 
tions 

The optical system used by the group at NPL is 
pictured in Fig. L Just to the left of plane I there is 
an integrating sphere of radius 25 cm which has 
a 450 W tungsten halogen lamp at its center, Plane 
I contains the exit aperture, <r, of the integrating 
sphere, a is a circular aperture of radius 
«o=0.12 cm, and it is in the front focal plane of a 
lens, L|, whose focal length is/i = 5 cm and radius is 
di=0.45 cm. 



Plane I 



Plane U Plane DI 
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The exit aperture of the integrating sphere is im- 
aged onto plane II. It is a straightforward matter to 
show that the magnification, M , for this imaging is: 



Af=-/afi=-4, 



(1) 



Two different kinds of things were done as regards 
plane IL In one case an aperture A of radius ai was 
placed in plane II, and experiments were done for 
the values a 3 =0.012 cm and 0.50 cm. 1 In the second 
case no aperture was used. For notational conve- 
nience, we will treat the latter case as if an aperture 
A of radius a$= « were in plane II. Also, since no 
shift was observed in the 0.50 cm aperture case, we 
will not consider that case* Plane III, where the 
measurements were made, is a distance z = 100 cm 
from plane IL 

The following notation will be used throughout 
this paper. The locations with respect to the optical 
axis of points in the planes I, II, and III will be 
specified by, respectively, the two-dimensional posi- 
tion vectors p", p\ and p. The locations of points in 
the planes occupied by Lj and Lj will be specified 
by, respectively, the two-dimensional position vec- 
tors ori and a 2 < The lengths of vectors will be de- 
noted by the corresponding non-boldfaced symbols, 
e.g., p w -\p% 

The following assumptions will be made. As con- 
cerns the properties of the light coming from the in- 
tegrating sphere, it will be assumed that the exit 
aperture of the integrating sphere radiates as a uni- 
form lambertian source* Let the spectral radiance 
of this source, at angular frequency m, be denoted 
byiJ^w). The cross-spectral density, W il) (pT, pj, o») 
of the light in plane I can be written as [9-1 1 J: 2 

WM(p'!,f&,*>) = 2<irB a (a>) 



xJoikltf-plftPfaDPoiri), 



(2) 



Fi& 1. The optical system. 

There is an interference filter, IF, a distance 
D\ =20 cm behind Li, and a second lens, L^ a dis- 
tance Di = 30 cm behind the interference filter. L 2 is 
an achromatic doublet, and its focal length is/ 2 = 2Q 
cm. The radius of La and the transverse dimensions 
of the interference filter are all larger than a x . It is 
therefore a straightforward matter to show that Lj 
is the aperture stop for the system. In their experi- 
ments the group at NPL used six different interfer- 
ence filters; the shortest peak transmission 
wavelength used was 422 nm and the longest was 
652 nm. 



where k is the wavenumber of the light, 

fc= < o{c = 27r}\ t (3) 

/o is the spherical Bessel function of the first kind of 
order zero. 



i In Ref. [1] the diameter of the smaller aperture was incorrectly 
reported as 0.24 cm. The correct value is 0.024 cm [8]. 
2 In order for Eq. (2) to be appropriate, the radius of lhe source 
must be much larger than the wavelength of the light In the NPL 
experiments the wavelengths of interest are visible and the 
radius of the exit aperture of the integrating sphere is 0.12 cm, 
so this condition is fulfilled, 
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;"o(m) = 



sin(u) 



(4) 



and Po is the pupil function for the exit aperture of 
the integrating sphere, 



Po(p")=ciic(p"/a ). 



(5) 



In Eq. (3) c is the speed of light in vacuo and A is 
the wavelength of the light. In Eq. (5) circ is the cir- 
cle function, 



circ(tf) = l, if |u|£l, 
=0, if M>1. 



(6) 



function we must first discuss the effect of the inter- 
ference filter. 

The effect of the interference filter on a plane 
wave is depicted in Fig. 2. Here a unit amplitude, 
monochromatic plane wave of angular frequency qj 
propagating in the direction specified by the unit 
vector s is incident upon the filter at an angle of in- 
cidence 0. If, as we have assumed, the finiteness of 
the transverse size of the filter can be neglected, the 
effect of the interference Filter is to change the am- 
plitude of the plane wave by the factor t(8 f a>) t the 
amplitude transmission function of the interference 
Filter. 



Secondly, since Li is the aperture stop for the sys- 
tem, it will be assumed that the finiteness of the 
transverse sizes of L-2 and the interference filter can 
be neglected. Thirdly, it will be assumed that the 
transmittance of the interference filter may be de- 
scribed by the Lissberger-Wilcock model [12-13] T 
and that its behavior is such that the approximation 
described in Sec. 4.1 is valid. In addition, both 
lenses will be treated as thin lenses, and the possi- 
bility of chromatic aberration introduced by either 
Li or L^ will not be taken into account. The paraxial 
approximation will be used throughout our work. 

3. Cross-Spectral Density of the Light in 
Plane II 

In this section* the cross-spectral density in plane 
II will be investigated. First the coherent impulse 
response function for the propagation from plane I 
to plane II will be determined. This result will then 
be used to obtain an "exact" (within the paraxial 
approximation) expression for the cross-spectral 
density incident upon plane II. Approximations ap- 
propriate to the NPL experiments will then be used 
to obtain an approximate expression for the cross- 
spectral density in plane IL The corresponding 
spectrum and complex degree of spectral coherence 
in plane II will then be calculated and discussed. 

3*1 Coherent Impulse Response Function for the 
Propagation from Plane I to Plane n 

Let P* be a point in plane I located by position 
vector p", and let P' be a point in plane II located 
by position vector p\ The coherent impulse re- 
sponse function, h{p\ p", <y), for the propagation 
from plane I to plane II is the field at P' due to a 
monochromatic, unit amplitude point source of an- 
gular frequency & at P*. In order to calculate this 



~w 



exp(iks.r) 



t(0,o>) exp(iks-r) 



IF 



Fig. 2. The effect of the interference filter on a plane wave. 

When this effect is taken into account, it can be 
shown (see Appendix A) that in the paraxial 
approximation h(p\ p", a>) is the product of t(p 7/:, 
w), the amplitude transmission function of the in- 
terference filter evaluated at 0=p7/2, and hx(p\ 
p", w), the coherent impulse response function 
(paraxial form) for the case when there is no inter- 
ference filter present, i.e., 



h{p\ p", a>)=t(p'ff 2t a>)h N {p\ p\ a>), (7) 



where 



Mp\ P ", ») = C( W )exp{i*W(p-) + x(p')]} 



xBesinc (^ | p"-^p' |) (8) 



and 



<*■>-»$ <*> 



*(p*)-/i + 0>" 2 /2/i). 0<» 
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Be S inc(«)=2/ 1 („)/, - (12) E «• * ->-«l{*(l-g )^-^)a fi ] , (16) 



In Eq. (11) D is the distance from Li to L2, i.e., 
D -Di +£> 2 . In Eq. (12) J v is the Bessel function of 
the first kind of order one. 

3.2 Cross-Spectral Density of the Light Incident 
Upon Plane II 

The cross-spectral density of the light incident 
upon plane II, JK (II_) (pi, #2, <o), is related to the 
cross-spectral density in the plane I by the expres- 
sion [14] 

W^-\p[, ri,«>)=j J k*(p' h p" u v)h{& p2, *>) 



xwM{p" u A w)d 2 pid 2 p?, 



(13) 



where h (p' t p", <u) is the coherent impulse response 
function for the propagation from plane I to plane 
II. Upon substituting Eqs. (2) and (7) into Eq. (13), 
and using Eq. (8) in the result, we find that 



W^->ipl, P i, 0>) = 2TTB^)t*(p' i /f 1 , «)f<pi// 2 ,a>) 



xF(p[,p{, w), 



(14) 



where 



Fip'upi, <o) = C(o))E(p\, p\, o>) 
xj $E (pi, A w) Besinc (^ | P I - ± p[\) 

a a 

xBesinc(^|p5 - ^ pi | );« (*|p8-rf|) 



E(pl /*, »)-exp[*(/« 2 -rf 2 )/2/.]. (17) 

Equation (14), with F(p\, pi, <w) given by Eq. (15), 
is an exact (within the paraxial approximation) 
expression for the cross-spectral density incident 
upon plane II. However, the integral on the right 
hand side of Eq. (15) cannot be evaluated analyti- 
cally. Nevertheless, for the NPL system certain 
approximations are appropriate which simplify 
F(p[, p{ t a>) considerably. 

First, it is shown in Appendix B that for the NPL 
parameter values we may set £(pi, pS, (o)^l in the 
integrand on the right hand side of Eq, (15), Hence, 
to a good approximation 



F(pU pi, <0) = C\n)E( P [, & 6>) J J 

tJ IT 

x Besinc(^| P T-^pi | ) Besinc (^i|rf - 1 pi | ) 



x/o(*|p5-pi-|)dVd 2 p5 • 



(18) 



Next, it is shown in Appendix C that for the NPL 
parameter values, 



1/ <*»»< (7rk-iM) 

fT IT 

x Besinc (^| pS-^/* | )jo (* |p5 -pl| )d*rfd 2 pS 



^""-MfrM) 



ic 



xd 2 p?d 2 p?, 



(15) 



x circ(p[/|M |ao) circ(pS/|Af \a„) . (19) 
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It then follows from Eqs. (14), (18), and (19) that 
the cross-spectral density of the light incident upon 
plane II is given, to a good approximation, by the 
expression 



77Cl 



xtitiffi, <o)E(pU pi, *i)Besino(^ |pi~pi ) 
xcirc(pI/|A/|fl )circ(pi/|M|a ). (20) 



The spectrum, 5 (It_) (p\ w), of the light incident 
upon plane II at the point P r located by the vector 
p' and at the angular frequency <o is obtained from 
the corresponding cross-spectral density by setting 
pi=P2=p'- J It therefore follows from Eq. (20) that 



j^B^Mp'fc, *,)Pcirc(p7|A/M. (21) 



xt 



M^W. Pi, «)= ^5 ( <u )f*(p 1 7/ 2j w ) 



(piffi, a>) E(p[, pi, o>) Besinc^y 1 pi -pi J 



xcirc(p{/a) drc(pifa), 



(23) 



where a is the smaller of |Af|a > and fl 3 , i.e., 
a = min()A/|ao T a^). 

The physical significance of the parameter a is 
that it is the radius of the "secondary source." If 
fla < \M\ao part of the incident light is blocked by the 
aperture A, and therefore the radius of the sec- 
ondary source is a± m However, if a$ >\M\ao the aper- 
ture does not (within the approximation we are 
using) block any of the light, therefore |A/|ao is the 
radius of the secondary source. The distinction be* 
tween these two cases is necessary because when 
the small aperture is used a* < |A/|ao, whereas when 
no aperture is used d 3 >lA/|flt>. 

The spectrum, S m (p*, w), of the light in plane II 
can be calculated by using the result in Eq. (23): 



It can be seen from this expression that this spec- 
trum is nonuniform (since / depends on p') and oc- 
cupies the region p'^jAfla^ i.e., the region which 
corresponds to the geometrical optics image of o-, 
the exit aperture of the integrating sphere. 

33 Cross-Spectral Density, Specrrum, and Com* 
plex Degree of Spectral Coherence in Plane II 

The cross-spectra! density, W w {p\, p5, w) in 
plane II (with the effect of the aperture A taken 
into account) is given by the expression 

IP™ (pi, /*,") = 
W ^ (Pi, p5» «) c\rc(plIai)tiTc(pya 3 ). (22) 



^(P'^) = ^" , (P , .P» ) 



7711 



- j£ JM*>)l'(p7/ 2 , o>)fcirc(p7a), (24) 



We see from this expression that the secondary 
source occupies a circular domain of radius a t and 
that the spectrum of the light inside the secondary 
source is proportional to the product of spectral ra- 
diance in plane J, /?u(ft>), and transmittance of the 
interference filter, \t(p'tf 2% a>)\ 2 . 

The complex degree of spectral coherence, 
P iU) (pu p5, w), of the light in the plane II is defined 
by the formula [16]: 



Upon substituting the approximate Eq, (20) into 
Eq. (22), we find that the corresponding approxi- 
mate expression for the cross-spectral density of the 
light in plane II is 



^(pi, ,*,*,) = 



0*">(pl. Pi, <») 



VS m QA, *>) VS m fa vY 



(25) 



3 The spectrum S h also sometimes referred to as the spectral 
density of the light. It h the optical intensity per unit angular 
frequency ai the position and angular frequency in question [15]. 



Upon using Eqs. (23) and (24) in Eq. (25), we find 
that when p[<a and fn <a, i.e., for pairs of points 
located inside the secondary source, 
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xBesincl 



($ w-4 



(26) 



This result shows that the group at NPL was incorrect 
in stating that the filter-lens combination eliminates 
the wavelength dependence of the complex degree of 
spectral coherence of the light in plane IL In particu- 
lar, since the first three factors on the right hand 
side of Eq. (26) are each unimodular, it follows that 
the modulus of ji tU) (pi, pj, o>) is equal to the abso- 
lute value of Bcs\nc(ka t \pi- p[\ff2)'. Hence, the as- 
sertion by the group at NPL that the wavelength 
independence of p {U \p\ t f&, w) causes the spectrum 
to change upon propagation from plane 11 to plane 
III is not correct. There may be a change of the 
spectrum (we will investigate that in the next sec- 
tion), but it is not caused by a wavelength indepen- 
dence of p iU) (pU pi, <a). 

In the next section we will refer to the effective 
correlation length, L (to), of the light inside the sec- 
ondary source. The correlation length of a Bestnc 
type correlation function is typically taken as the 
smallest separation, (p2-p[[, for which the complex 
degree of spectral coherence takes on the value 
zero. Since Besinc has its first zero when its argu- 
ment is equal to 3.832, it follows from Eqs. (26) and 
(3) that 



L(a})= 0.610 






(27) 



4. Spectrum of the Light On-Axis in 
Plane III 

Let P be a point in plane III located by position 
vector p. When the Fresnel approximation is used 
to propagate the field from plane II to plane III, the 
spectrum of the light at P is given by the expression 
[17] 

S™(p t m) = 

yff \l ^ Cn) (Pi ti> *>) expPW~ P T 2 )/2z] 

xexp[-^(pj-pl) •p/*]dVid a pi, (28) 

where W m (pU pS, &>) is given by Eq. (23), and the 



integrations extend formally from - » to °o in each 
variable. 4 

If the secondary source obeys the Leader condi- 
tion [18] 



kaL(oi) 



«1, 



(29) 



then we are in the far zone of the secondary source 
and the quadratic exponential term in Eq. (28) can 
be dropped. We will restrict our considerations to 
the on-axis observation point, p =0. 

4.1 No Aperture in Plane II 

In this case a — ]A/|ao, and the Leader condition is 
not fulfilled for the NPL parameter values. The on- 
axis spectrum is therefore given by the expression 

5<m,(0 ' a)= jkf- W w ' <l,)(p1, ** <w) 

x exp[iJfc (tf 2 - pf*)f2z}&pi d 2 p5, (30) 



where W m (pi, pj, to) is given by Eq. (23), with 
a = |A/|ao. 

In order to investigate ^'"'(O, &>), let us first in- 
troduce average and difference vectors, p' and />', in 
plane II: 



p'=Hri + 0l), 



P'=ti-Pi 



(31a) 
(31b) 



In terms of these variables, it follows from Eqs, 
(30), (23), (16) and (27) that 

S m (0, W )=- x ^ 2 5o(«)|J G(p\ p\ a>) 

x exp( -ikp' • p'/d)d 2 p 'd 2 p', (32) 

where 

G(Jp',p', <0)=t*(\p'- l 2 p'\/f2, o>) 

xrflp' + kp\/fi, *>) Besinc[3.832p7L(o>)] 
Xcirc[(p'-|p')/|A/lfl„]circ[(^'+ip')/|A/|a fl ] , (33) 



4 In practice, lhe circ functions in Eq, (23) limit the domains of 
integration to the area occupied by the secondary source. 



272 



Volume 99, Number 3, May-June 1994 

Journal of Research of the National Institute of Standards and Technology 



and 



d h I/: 1 ) z 



(34) 



Let us now examine the behavior of G(p\ p, a>) 
as a function of the difference vector p' for a fixed 
value of the average vector p'. The region over 
which the Besinc function is significant is the do- 
main p'£L(a>). As p* explores this region the first 
arguments of r* and t in Eq. (33) each change by at 
most the angle 



M^ = 0.610 A 



(35) 



For visible light this is an extremely small angle. For 
example, for A =550 nm and ^i = 0,45 cm, L(o>)/ 
fi-15 uxad. Therefore we will assume that the be- 
havior of the interference filter is such that when 
evaluating the integral on the right hand side of Eq. 
(32) we may, to a good approximation, neglect the 
p' dependence of both t* and r, i.e., that we may 
approximate G as 

G(p\ p\ a>) = |f (pV/i. ")l 2 Besinc[3.832p7L(«)] 



xcirc[(p f -ip')/|A/|flo]circ[(p' + ^p')/|A/|anl (36) 



Furthermore, for the NPL system parameter values 
and the wavelengths used in their experiments, 
L(to)«\M\a . Therefore we will now use the 
quasihomogeneous approximation [19] which corre- 
sponds to replacing Eq. (36) by the expression 

G(p',p\o>) = 

H~P*lfa 6>)] 2 Besinc[3.832p'/L(w)] circ(p7|A/|flo). 

(37) 

Upon substituting the approximate form Eq. (37) 
for G into Eq. (32), it can be shown, after some 
straightforward calculations, that the on-axis spec- 
trum in plane III is given by the expression 

S m (Q> w)=^Bo(co)M(<o) , (38) 

where 



A/(^)=^ J \t$'ff Zr 6>)p cire(p7|M|flo) 



xcirc(p7/0dV , 



(39) 



and b —a\dlfi. For the NPL parameter values* 
|iVf|tfo^0.48 cm and 6=0.35 cm. Therefore, Eq. 
(39) simplifies to 



M(<*)= ^i J|r(P'#, o>)\ 2 circ(p7A)d 2 P\ 

2 r 
= j 2 jjt(e.a>)\ 2 eae, (40) 



where $=blfi. 

Let us now discuss the transmittance, |/(0, tu)\% 
of the interference filter. We will assume that it can 
be described by the Lissberger-Wilcock model [12- 
13], This model works well for Fabry-Perot and all- 
dielectric interference filters for angles of incidence 
which are less than 20°. It can also be used to de- 
scribe, less accurately, double halfwave and induced 
transmission filters [20]. 

Let Ao denote the peak wavelength transmitted by 
the interference filter at normal incidence (0=O) t 
(AA) D denote its bandwidth (FWHM) at normal in- 
cidence, and 7o denote its maximum transmittance. 
According to the Lissberger-Wilcock model the 
transmittance of the interference filter is given by 
the expression 



I'f'-'M'+P^ + ^n 



-1 



(41) 

where <n is the effective index of refraction of the 
interference filter. 

The integral in Eq. (40) can be evaluated (see 
Ref. [12]) without any approximation, for the trans- 
mittance function Eq. (41). The result is that 



"(^'"-'[ i+fg+ao ]' < 42 > 



where X = Ao/? 2 /[2tt(A A >,], and £ = 2(A - Ao)/(AA) . 
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Some comments about the interference filter 
parameters Ao, (AA)o, and t\ are now in order. The 
effective indices of refraction, i.e., the values of tj, 
for the filters used in the NPL experiments are not 
known. However, typical values of r\ are in the 
range 1.4<ti<3.4, with the lower values being 
more typical [21-25], We will therefore use v\ = 2 in 
our calculations. The precise values of Ao and (A A)o 
are also not known. For reasons which will be ex- 
plained in Sec, 4,3, we will take the values of Ao for 
the six filters to be 421.9 nm, 484.0 nm, 512.4 nm, 
566,0 nm, 609,1 nm, and 652.0 nm, with their band- 
widths being, respectively, 9 nm, 9 nm, 5 nm, 13 nm, 
8 nm, and 8 nm. 

4*2 Small Aperture in Plane II 

In this case a$- 0.012 cm and the Leader condi- 
tion is fulfilled. It therefore follows from Eqs. (28) 
and (23) that the spectrum can be written as 



£< m >(0,<o) = 



vt*l 



TUi 



kh*f 2 



Kl B^m) J J t*(p[/f 2 , m)t(pi!fi, to) 

A A 



xE(pl, & <*) Besinc(^ |p£-pi|)d 2 pid 2 pL (43) 



where A is the circular aperture of radius a^ For the 
same reasons as were discussed in the Sec, 4.1, we 
will assume that the product of r* and t on the right 
hand side of Eq. (43) can be replaced by |/| 2 evalu- 
ated at the average position, i.e,, that to a good ap- 
proximation, 

£f ]]I >(0, «)- 






fl (») J J K|p! + pi\/2f it o>)\ 2 E(p\, pi, to) 



x Besinc 



(t?w-*i) 



dVidVS, 



(44) 



For the NPL parameter values, it can be shown 
from Eq. (41) that, to a very good approximation, 
|'(|p5 + pJ|/2/2, *>)IH<(0, *)\\ for all p| and & 
which are in A. The on-axis spectrum in plane III is 
therefore given* to a good approximation, by the ex- 
pression 



S<i>i> (0j w)== ^3 5o(w)]f(0j w) p jv(„) f (45) 



where 



"M-I& J /*<**•> 



xBesincfe |pJ-pJ|) d 2 pld 2 ^ . (46) 

This integral can be evaluated by substituting the 
Fourier integral representation of the Besinc func- 
tion [see Eq. (88)] into the integral, interchanging 
the orders of integration, and then recognizing the 
resultant integrals as familiar diffraction integrals. 
The result is that 



*(*>;= i-i^^fa*^), (47) 



where 



kaUD A 



v = - 



IT' 



(48) 



(49) 



and the Qi%{v) functions are those introduced by 
Wolf [26] and simplified by Petersen [27]: 



ft(v)^(v)+;i(v), (so) 



<M») = - 2^71) {v[y2(v)y2s+,(v) 



-ii(vy i+1 (v)]+2sy,(vy li+ i(v)} f if j>o. (si) 



A3 Numerical Investigation of the Spectrum On- 
Axis in Plane III; Comparison to the NPL 
Experimental Values 

Let s$ IJ} (0, A) be the spectrum, as a function of 
wavelength, at the on-axis position in plane 111, for 
the case in which there is no aperture in plane 1L 
Let the wavelength at which this spectrum peaks be 
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denoted by A P and the bandwidth (FWHM) of this 
spectrum be denoted by AA. Let 5ji ,l> (0, A) be the 
spectrum, as a function of wavelength, at the on- 
axis position in plane III, for the case where the 
small aperture is in plane IL Let A P denote the peak 
wavelength for this spectrum. The group at NPL 
measured Ap, AA, and A P for six different interfer- 
ence filters [1]. Their results are shown in Table 1, 
along with the shift, 

5A p =Ap-Ap, (52) 

which occurred due to the insertion of the aperture. 

Table 1. NPL experimental results 



Ap, in am 


AA 


, in nm 


Ap, in nm 


SAk 


, in nm 


422.0 




9 


421.0 




-1.0 


484.1 




9 


483.6 




-05 


512.4 




5 


514.1 




1,7 


566 J 




13 


564.1 




-2.0 


609.1 




8 


610.3 




1.2 


652.0 




8 


653.2 




1.2 



There are several things which should be noticed 
about the results shown in Table 1. First, the shifts 
are* in absolute value, of the order of 0.5 nm to 
2.0 nm. Secondly, some of the shifts are blueshifts 
and some are redshifts. Furthermore, there are no 
obvious trends in the behaviors of the shifts as func- 
tions of either A P or AA. 

Let us now investigate the shifts that our theory 
predicts. It follows from Eqs. (38) and (42) that the 
spectrum when there is no aperture in plane II is 
given by the expression 

*n0, A) = 4\(A)m N (A), (53) 



where b^k) is the spectral radiance, as a function of 
wavelength, at the exit aperture of the integrating 
sphere and 



ira\ 



*r(0,A)=^& (A>n A (A), 



(55) 



where 



m A (X) = T N(2ijc/\)\\ + 



/A){l 



L (AA)o J I 



(56) 



In order to simulate the experimental conditions, 
we did the following things. The color temperature 
of the lamp was 3200 K [7], We therefore took the 
spectral radiance, jE>o(A), to be a Planck distribution 

[28]: 



VA) = 



2hc 2 



1 



A 5 exp(hc/XknT)-l 



(57) 



with a temperature of 7=3200 K. In Eq. (57), h is 
Planck's constant and &b is Boltzmann's constant. 
Furthermore, we chose the center frequencies, Ao, 
and the bandwidths, (AA)o , of the interference fil- 
ters to be such that the Ap and AA values obtained 
for the spectrum ^Jg Tr) (0, A) [from Eq + (53)] agreed 
with the experimental values of Table 1 to the num- 
ber of decimal places being displayed there. 

Table 2 lists the interference filter parameters, 
and Ao and (AA) U , we used in our calculations and 
gives the shifts, SA Pj predicted by our theory. These 
shifts were obtained by using a search routine to 
find Ap and Ap from, respectively, Eqs. (55) and (53), 
and then subtracting the two values. 

Table 2, Theoretical values for the shift 

Ad* in nm (AA)n, in nm &Ap* in nm 



421.9 
484.0 
512.4 
566.0 
609. L 
652.0 



9 
9 

5 
13 

8 
8 



0.007 
G.0O9 
0.009 
0.010 
0.011 
0.012 



^'"I iot I' 



(54) 



Here f and X are given by the expressions which fol- 
low Eq. (42). It follows from Eqs. (45), (47), and 
(41) that the spectrum when the small aperture is in 
plane II is given by the expression 



These shifts disagree with those measured by the 
group at NPL. The theoretical shifts are approxi- 
mately two orders of magnitude smaller than the ex- 
perimental shifts. Furthermore, the theoretical 
shifts are always redshifts, whereas the experimen- 
tal shifts are in some cases redshifts and in other 
cases blueshifts, 
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At this point some further comments about the 
shifts predicted by our theory are in order. First, let 
us consider Ap* £ ^(O, A) is the product of two wave- 
length dependent factors, m A (A) and 6 (A). Due to 
the presence of the factor N(2ttc/X), m A (A) is not 
centered at Ao, it is blueshifted from it by an amount 
of the order of 0.0001 nm. Multiplication by the 
Planck spectrum then causes m A (A) to be redshifted 
by about 0.1 nm. Now, let us consider Ap. sft u) (0, A) 
is the product of two wavelength dependent factors, 
Wn(A) and io(A). win(A) is not centered at fa, it is 
blueshifted from this value by an amount of the or- 
der of 0-01 nm. Multiplication by the Planck spec- 
trum then causes jhn(A) to be redshifted by about 
0.1 nm. When we subtract Xp from Ap the two Planck 
shifts cancel. Since the blueshift of mti(\) is much 
larger than the blueshift of m A (A), the resulting 
shift, §Ap, is a redshift of the order of 0.01 nm. 

The physical origin of this redshift is as follows. 
The field incident upon the interference filter can 
be represented as a superposition of polychromatic 
plane waves traveling in different directions, and 
the wavelength at which the transmittance of the in- 
terference filter peaks decreases as the angle of in- 
cidence of a plane wave increases [see Eq. (41)]* 
Therefore, for each polychromatic plane wave inci- 
dent upon the interference filter at an oblique angle 
of incidence, the spectrum of the transmitted light 
has a peak wavelength which is less than it would be 
for a normally incident plane wave. In the case in 
which no aperture was used, the contributions from 
a significant set of such waves arrive at the on-axis 
observation point in plane III, However, in the case 
in which the small aperture was used, the contribu- 
tions from a large subset of these waves were 
blocked at plane II and did not arrive at plane III. 
As a result, the spectrum in plane III when the 
small aperture is used peaks at a longer wavelength 
than it does when no aperture is used* 



5. Conclusion 

Our paper contains two separate sets of results. 
(1) In Sec. 3 an approximate form for the cross- 
spectral density of the light in plane II was obtained 
[Eq. (23)], and it was shown that the corresponding 
complex degree of spectral coherence contradicts 
the explanation for the shifts given in Ref« [I J (2) In 
Sec. 4 further approximations were made to propa- 
gate the cross-spectral density obtained in Sec. 3 
from plane II to plane III, and an approximate ex- 
pression for the on-axis spectrum was obtained, 
both for the case in which no aperture is used, and 
for the case in which the small aperture is used. It 



was found that the peak wavelength of the spectrum 
in the latter case is shifted with respect to the peak 
wavelength for the former case. However the shifts 
predicted by our analysis are much smaller than 
those reported in Ref* [1], so small as to be unob- 
servable to within the accuracy of their experi- 
ments. 

This brings us to an important question. Our 
analysis is predicting no observable shift (to within 
the accuracy of the measurements of the group at 
NPL), and yet the group at NPL observed shifts; so 
where is this shift coming from? In our opinion 
there are three possibilities. One possibility is that 
their interference filters did not behave in the man- 
ner we have assumed in our calculations, and that 
shifts arose as a result. A second possibility has to 
do with the spatial coherence properties of the light 
in plane III. Since the secondary source created in 
plane II when no aperture is used has a radius 
which is forty times as large as in the case in which 
the small aperture is used, the spatial coherence 
properties of the light in plane III will be quite dif- 
ferent in the two cases. If this difference is signifi- 
cant enough, it may be that the monochromator 
used in the detection process responds differently 
in the two cases. A third possibility, and in our opin- 
ion the best one, is that the shift is caused by chro- 
matic aberration introduced by Li, or L;, or both. Li 
has a short focal length and is not an achromatic 
doublet, and is therefore an obvious candidate for 
introducing significant chromatic aberration. L 2 is 
an achromatic doublet. However, it can only be per- 
fectly achromatic for two wavelengths, whereas fil- 
ters with six quite different center wavelengths were 
used in the NPL experiments. 

Finally, let us comment on the relevance of these 
possibilities for typical spectral irradiance measure- 
ment systems. Typical spectral irradiance measure- 
ment systems do not use interference filters. There- 
fore, if the shift is caused by a nonideal behavior of 
the interference filter, it is irrelevant for typical 
spectral irradiance measurements. As concerns the 
second possibility, it has been known for some time 
that the spatial coherence properties of the light in- 
cident upon a monochromator do effect its re- 
sponse, viz., its slit scattering function [29-30]. 
Hence any comparison of experimentally measured 
spectra which have significantly different spatial co- 
herence properties at the entrance aperture of the 
monochromator should take this into account. Fi- 
nally, if the shift is caused by chromatic aberration, 
it will not occur in spectral irradiance systems which 
use only mirrors, and is an effect which should obvi- 
ously be taken into account in any system which 
uses lenses. 
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6. Appendix A: Coherent Impulse Re- 
sponse Function for the Propagation 
from Plane I to Plane II 

6.1 Field Exiting Li 

Let Q be a point in the plane occupied by Li, and 
let a\ be the two-dimensional position vector in the 
plane which locates Q. The field, U il ~ } (a u o>), inci- 
dent upon Li at the point Q due to a monochro- 
matic, unit amplitude, point source at position P" in 
plane I is given by the expression 



^ 1 -^.")-n' BtP |nl^ Pl) - < 58 > 

where ri = (a if /i). Since the angles involved are 
small, this expression may be approximated 
(paraxial approximation) as 



£/<'->(«,, „ ) = e *P(y/0 exp[ifc(« l -/>' l ) a /yi]. (59) 

The field exiting the lens, U^^\a u co), is given by 
the expression [31] 

t/<'->(ai, < U )exp(-i*oi/2/ 1 )P.{a'i) , (60) 



where P(ai) is the pupil function for Lu 
Pi(ai) =circ(«i/fli). Using Eq. (59) in Eq. (60) yields 



£/ tl+J (a„*0 = 



eXP[/ fv. (P ' )] «*-'** ■'"*)*<■») ' < 61 > 



where 



*(P")=/1 + (P"V2/,) 



(62) 



In the next subsection we will need U {l *\K , <o), 
the two-dimensional spatial Fourier transform of 
t/< 1+ >(«., «), 



0"* } (K, «)- 
SU tu \a u w) exp(-/tf • o,)d 2 a,, 



evaluated at Jf =ksi. It follows from Eqs. (61) and 
(63) that 

M/l 



- ^ exp[ift*(p")] 



x Besinc[ifliJ5i-5o±|] , (64) 

where s§± = —p^lfx. 

6.2 Field Incident Upon La 

Let P be a point in-between Li and the interfer- 
ence filter. Let r be a position vector which goes 
from the center of L] to P. The field at P can be 
represented as a superposition (angular spectrum) 
of plane waves [32]: 

U(r t a))=f a($ X1 o))exp(iks *r)d 2 $±. (65) 
where s=(s ±J s z ) is a unit vector, 



s^Vl-sl , ifJi^l, 



and 



='V*i-l » if$i>l, 



a(5, lW )=pi/ (l+J (h 1)W ). 



(66) 



(67) 



(63) 



In Eq. (66), ja = |sx|. 

In spherical polar coordinates, as measured from 
the center of L u *-(i, S, <£). In this notation 
s± = (Sx, $,) = (sin0cos4>, sin0sin<£). Therefore the 
angle of incidence of the plane wave expfjks -r) 
with respect to the normal to the surface of the IF 
is [see Fig. 2], where 5j_ = sin &**6. 

Let t(9, o>) be the amplitude transmission func- 
tion of the Filter for a plane wave of angular fre- 
quency wand angle of incidence 0. If we neglect the 
finiteness of the transverse size of the IF, the effect 
of the interference filter is to change the amplitude 
of each plane wave, exp(iks •#*), by the factor t(su 
w). The amplitude, call it a(s Xi &) t of the plane 
wave exiting the filter is therefore 
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d($x t a>)= t(s u <i))a($i, o>) t 



1 



h(p\p» t <o) = 



= -pt(si, w) U il+) (ks Xt a) , 



TTOi 



X Besinc[ftfli|s±-5oj.|] . 



(68) 



Let ri—{a2 y D) be the position vector which goes 
from the center of Lt to the point in L2 located by 
the two-dimensional position vector 02, The field, 
£/ (2 "*(ak, <u) incident on L2 can be represented as a 
superposition of plane waves: 

U {2 ~ ) (a 2 , to)= J a(s u o))exp(iks *r)d 2 s ±i 

= / a(s ±t w) exp(ik$ ± • 02) 

x exp(iks t D)d 2 s±. (69) 

Since the angles involve d are small, we may make 

the approximations =Vl —Sx**^ "l-^i m Eq. (69), 
and we find that 

U {2 ~\eti, (o) = txp(ikD)fi($u c*>) cxp(iksi- ai) 

x exp( - ikDs\l2) d 2 s ± . (70) 



ira\ 



-p^exp{iA:W(p-)+^{p')J 



xBesmc(^a,|(p'y/2)-0>7/i)|)f(p7/2, a>) , (74) 
where 

X (p-) =/ 2 +Z> + ( 1 -| ) (p' 2 /2/ 2 ) . (75) 



The magnification of the imaging from plane I to 
plane II is M = —filfu It therefore follows from Eq. 
(74) that the coherent impulse response function 
can be written as 

h(p\ p\ &)=t{p'tf 2 , oj)h H (p\ p\ to), (76) 

where /zn(p'jP*i fr>) is the coherent impulse re- 
sponse function if there is no interference filter, 



h»(p\ p\ a>) = C(to)exp{ik[<t>(p") + x(p')'}} 



xBesmc(^|^-i^|). (77) 



6J3 Final Expression 

The coherent impulse response, h(p\ p" 7 m), for 
the propagation from plane I to plane II is the field 
in plane II due to the point source in plane I, i.e., 



h(p' t p a ,w) = V™{p\u). 



(71) 



If we neglect the finite transverse size of Lz, this 
equation can be written as [33] 

ft(gV^«)° e ^ (72) 

where 4f(p')=f 2 + (p' 2 f2f 2 ). It follows from Eq. (70) 
that 

U i2 -\ks u «)*= 

\ 2 a (s ± , <u) cxp(ikD) exp(~ikDs 2 ±/2). (73) 

Upon substituting Eq. (73) in Eq. (72) and Eq, 
(68) in the result we find that 



Here $(?") and x(p') are given by Eqs. (62) and 
(75) respectively, and 



7. Appendix B: Derivation of 

Let us first define average and difference vari- 
ables in plane I as, respectively, 



p'=pj-pr. 



(79a) 
(79b) 



In terms of these variables E(p", pi, <o) [see Eq. 
(17)] can be rewritten as 



£(pT,p5.«)-exp(^r £"■•*•) 



(80) 
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As pi and p? each explore the domain n-, the max- 
imum value of p" = |p*| is a . In the integrand of Eq. 
(15) the range of values of p"=f&— pi over which /<> 
is significant is p"£A/2. Hence as pi and p£ range 
over the domain where the integrand is significant, 



2lT ^ ^ 



TTUq 



/. 



(81) 



For the NPL system aa = CU2 cm and fi=5 cm; 
therefore over the domain where the integrand is 
significant the real and imaginary parts of £(p7, p£, 
«*») fulfill, respectively, the inequalities: 



0,997 £ Re[E(pr, p5, <»)] < 1 , (82a) 

-0.075 < Im[£(p:, pj, «)] £ 0.075 (82b) 

8* Appendix C: Derivation of Eq. (19) 
Let 

B(pU & <u) = J|Besinc(^p |rf-^pi|) 

a a 

(83) 
and let us make the change of variables 

*j=t)*fa*, ti=h2) (84a) 

iy=«7Afa. > , = 1.2) (84b) 

In terms of these new variables, Eq. (83) can be 
rewritten as 

£(pl.P2» o>)=do J I Besinc(ici|i*i-V||) 

x Besinc(Ki \u 2 -v^{k \u 2 - «j|)d 2 w, d 2 i/ 2 , (85) 

where uA is the unit disk, u.d, ={«;u<1}, and 

kaia 



*1 = 



/i ' 



(86a) 
(86b) 



If ki » 1 then the first and second Besinc func- 
tions in the integrand on the right hand side of Eq. 
(85) are, respectively, very sharply peaked about the 
values U\~v\—p\IMa\s and U2 = V2—pyMao. As a re- 
sult, if vj or v 2 lies outside the unit disk, then the 
integral on the right hand side of Eq. (85) is approx- 
imately zero, and if V] and v 2 both lie within the unit 
disk the integrals in Eq. (85) maybe approximated 
by extending the limits to ± », i.e., 



B{pu piy fl>)wfl! circ(pf/Affl ) circ(pi/Affl ) 
x/J Besinc(tfi|uj-i>i[) Besinc(Ki|w2— H) 

x;oM w 2-Ui|)d 2 w,d 2 W2 ♦ (87) 



The right hand side of Eq. (87) may be simplified 
by using the Fourier representations of the func- 
tions in the integrand: 

Besinc(ici|ii;-i>J)=— 5 J circ(J&/#rj) 



xexp{-iKj • fa -v;)] d 2 K h (j = l t 2) y (88) 



1 f tiic(JT/«) . te-uflfK 



V1-(K/k) 



(89) 



Upon substituting Eqs. (88) and (89) into Eq. (87), 
we find, after performing some simple integrations, 
that 

5 (pi, P2, a>)=D(a>) circ(pJ/A/fl ) circ (pi/Mao) 
x f circ(*/«,) circ(^K) ^g.^^ 

-D((a) circ(p|/Mfl ) c\rc{pHMaa) — ;■ **' v i 

J Vl-(^/f) 

xexpt/K-^-vO^/C, (90) 



where 



°H-(^?, Pi) 
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and the step performed in Eq. (90) follows from the 
fact that Ki/fc=tfi//i<l. 
If (ai/fi) 2 << 1, Eq. (90) can be approximated as 

B(p r u pi, cj)=D(a)) circ(pi/Moo) circ(p£/A/a ) 

x J circ(AY«: l )exp[7A' • (v 2 -vj)d 2 A:, 

-D(o)) circ(pi/A/flo) ciTc(p2fMad)TrKi 2 

xBesinc(Ki|v 2 -V||) ? 



xcirc(/»2/|A/|(io) 



(92) 



For the NPL system, flo=0.12 cm and arffi = Q.09, 
and the peak wavelengths of the interference 
Filters used in their experiments were in the range 
422 nm to 652 nm. For these parameter 



values, 1.04xl0 3 <; 



/t 



<1.61x10 3 and (a&y 



= 8,10x 10~ 3 ; therefore, jo » 1 and {a x lf\) 2 « 1. it 
then follows that the approximation Eq. (92) is ap- 
propriate for the NPL system. 
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